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We analytically calculate the energy spectrum of a circular graphene quantum dot with radius 
R subjected to a perpendicular magnetic field B by applying the infinite-mass boundary condition. 
We can retrieve well-known limits for the cases R, B oo and B 0. Our model is capable 
of capturing the essential details of recent experiments. Quantitative agreement between theory 
and experiment is limited due to the fact that a circular dot deviates from the actual experimental 
geometry, that disorder plays a significant role, and that interaction effects may be relevant. 
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I. INTRODUCTION 

Shortly after its discovery four years ago jlj , graphene 
triggered tremendous research both theoretically and ex- 
perimentally 0, The excitement was partly created 
by the fact that charge carriers in graphene are described 
by the Dirac equation for massless particles and a linear 
dispersion relation. In particular, there is no energy gap 
between valence and conduction band. Impressive ex- 
perimental results like the unconventional quantum Hall 
effect d, 0] and Klein tunneling Q were achieved and 
can be explained with these special properties. First ex- 
periments on graphene quantum dots were carried out to 
study their energy spectra in a perpendicular magnetic 
field 0, Q . The understanding of their spectra, especially 
how they may differ from a conventional 2D electron sys- 
tem described by the Schrodinger equation, is emerging 

Here we look at the energyspectrum obtained by mag- 
netic field spectroscopy in [8| from the theoretical side. 
We derive an analytic expression for the energy spectrum 
in a perpendicular magnetic field in the first part, look at 
different limits in the second part, and compare theory 
and experiment in the third part of Sec. II. We conclude 
with some final remarks in Sec. HI. 



II. ENERGY SPECTRUM IN A 
PERPENDICULAR MAGNETIC FIELD 

Derivation 

We start from the free Dirac equation expressed in 
cylindrical coordinates and include a magnetic field ori- 
ented normal to the graphene sheet. We use the symmet- 
ric gauge for the vector potential, A = B/2{—y,x,0) = 
B/2(— r sin0, r cos 0, 0) with (j> being the polar angle. 
The Hamiltonian then reads 



and the Dirac equation is HiJ;(r, (p) = Eip{r, (j)) with the 
wave function being a two-component spinor, 'il){r,(f)) = 
("01 (r, (f)), (j))). The charge of an electron is given by 
— e, vf is the Fermi velocity, and cr = (ax, cry) are Pauli's 
spin matrices in the basis of the two sublattices of A- and 
B-atoms. The electron spin is neglected in our analysis. 
A mass-related potential energy V{r) is coupled to the 
Hamiltonian via the CTz Pauli matrix. The mass in the 
dot is zero, V{r) = for r < R, but tends to infinity at 
the edge of the dot, V{R) — + oo. In this way, charge car- 
riers are confined inside the quantum dot. This leads to 
the infinite-mass boundary which yields the simple con- 
dition that /V^i = Tiexp[i0] for circular confinement 
flU . We are aware that the boundary condition can 
be different in the experimental situation. The parame- 
ter T takes the two values ±1 which leads to a distinction 
between the two valleys K and K' described by the Dirac 
equations of Eq. ([T]). Hence, in the following we can set 
V{r) — and we will respect the different energy spectra 
of the K- and K'-valleys via the boundary condition. 

Since the operator for the total angular momentum, 
Jz = Iz + ^cFz commutes with H, [H,Jz] = 0, we can 
construct simultaneous eigenspinors for H and Jz (m be- 
ing an integer). 
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Plugging this expression into the Dirac equation and de- 
coupling the system of differential equations, we arrive 
at a second-order differential equation for, e. g., Xii''') 
which depends only on r. 
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Xi{r)=Q. (3) 



H = VF (p + eA) ■ a + TV{r)(Jz, 
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The energy E is related to the wave vector k according to 
E = hvFk. We have introduced the magnetic length Ib — 
^ fij (eB). In order to solve this differential equation, we 
make the ansatz = exp [— r^/4Z^] ^ (r^)- This 
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yields the associated Laguerre differential equation 
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e (^1 = 0. 
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with f := r^. The solution is ^ (f) = 
cL (/c2;|/2 — (to + 1), m, f/2/|) , where L{a,b,x) is 
the generalized Laguerre polynomial and c is a nor- 
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malization constant. The second linearly independent 
solution of Eq. (|4|), the confluent hypergeometric 
function of the second kind, does not appear in the wave 
function ^ because it cannot be normalized. With the 
final result for ^ and hence for -01, -02 can be derived 
directly from Eq. ((T|). The wave functions finally read 



Mr. <t>) = cie^("+i)"^r"e-'''/'*'' 



2;2 



kH 



2 

r/lB 



kU 



{m + l),m,r^/2l% 



L 



(TO + 2),m+l,rV2?|) + L' ^'^^ 



(to+ 1),to, r^/2/ 



(5) 



Employing the boundary condition discussed above, we finally get an implicit equation for determining k, namely 
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Since the generalized Laguerre polynomials are oscilla- 
tory functions, there is an infinite number of knS for 
given B, to, and r which fulfill the above equation. This 
defines the radial quantum number n which labels the 
roots of the left part of Eq. ©. Here, we restrict our- 
selves to positive solutions for k so that ki > 0. There- 
fore, the energy spectrum E{n, m, r) of electrons confined 
to a circular graphene quantum dot which is exposed to a 
perpendicular magnetic field is determined through Eq. 

The relation —E{n, m, r) = E{n, to, — t) is a mani- 
festation of the electron-hole symmetry. 



Limits for B — > and R/Ib —> oo 

For a better understanding of Eq. ([5]), we will look at 
the two limits B ^ and R/Ib oo separately. Bessel 
functions of the first kind can be expressed as a limit of 
the generalized Laguerre polynomial , 



lim 

a — >oo 



1l (a, 6,^)1 = x-^/Vb (2VS) 



Using this property, Eq. ^ can be simplified to 

Tj„r{kR) ^ Jm+l{kR) 



(7) 
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in the limit B 0. This is the result already derived 
in The relation can be used to estimate the num- 



ber of charge carriers confined on a graphene dot when 
the energy of an excited state is measured 0, @| . More- 
over, we can deduce that E{n, to, r) — E(n, —to — 1, — r) 
for i? = 0. This is derived from the property Jm{x) = 
(— 1)™ J_m(a;). Hence, pairs of states are degenerate at 



zero magnetic field. There is no state at zero magnetic 
field and zero energy llj. This leads to an energy gap 
between negative and positive energy states. The size of 
the energy gap will be discussed below. 

Landau levels should be retrieved from Eq. ([6]) if the 
confinement is lifted. Mathematically, this is achieved for 
R/Ib — > oo. We express the generalized Laguerre poly- 
nomial in terms of the confluent hypergeometric function 
of the first kind M(a, /3, 7) [H 



L(a, b, x) 



M{-a,b+ l,x). 
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For R/Ib 00 , a power expansion of AI is possible and 
yields to first order 
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(i + o(i7r^)). (10) 



Rewriting the binomial coefficients with Gamma func- 
tions r(a;) and using one of their defining relations, 
r(x -I- 1) — xV{x), algebraic manipulations of Eq. ^ 
give 



E,ri — hvpkm = ivpy^2ehB{m + 1). 



(11) 



Hence, we retrieve the well-known Landau levels for 
graphene. Therefore, there will be a transition, governed 
by the parameter R/Ib, from a regime where the energies 
of the electrons are dominated by confinement (Eq. ^) 
to Landau levels (Eq. (fTTj) ). This transition including 
these two limiting cases is described by Eq. 
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FIG. 1: (Color online) Energy spectrum of a quantum dot 
with R = 70 nm. The formation of the lowest Landau levels 
can be seen as predicted by Eq. (|lip . Energy states for 
r = +1 are drawn with black, solid lines, those for r — — 1 
with green, dashed lines. 



Numerical Results and Comparison to Experiment 

We evaluate Eq. ^ for a dot of radius R = 70 nm 
which is about the same size as the device measured in 
Q . The energy spectrum as a function of magnetic field is 
shown in Fig. [T] for m — —4, . . . , 4 and n = 1, . . . , 6. For 
B = 0, the energy states are not equidistant. For higher 
magnetic fields, we can see the formation of Landau levels 
as it is expected according to the previous discussion. 
The zero energy Landau level is formed by states with 
quantum number t = — 1 and E > and those with 
T — +1 and E < 0. For completeness, we plot the first 
negative energy states as well in Fig. [TJ 

The lowest positive energy state has an energy of about 
4meV for B = 0. This gives an energy gap of around 
8meV between electron and hole states. Since the en- 
ergy gap to the next excited state is much lower, the 
electron-hole transition may be detected experimentally 
by a confinement enhanced energy. 

For quantum dots in semiconductors, the Darwin-Fock 
model is often used to qualitatively explain the experi- 
mental observations. In contrast to the model presented 
in this paper, the Darwin-Fock model is based on har- 
monic confinement giving rise to equidistant and highly 
degenerate energy levels at _B = 0. However, the general 
evolution from single-particle states at S = to Lan- 
dau levels at high magnetic fields is similar in the two 
scenarios. 

We recently performed transport spectroscopy mea- 
surements on Coulomb blockade resonances in a graphene 
dot with a radius of about 70 nm Fig. [2] is the main 
result of that experiment. It shows the position of con- 
ductance resonances as a function of magnetic field. The 
vertical energy axis was obtained by converting plunger 
gate voltage into energy using the measured lever arm. 
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FIG. 2: (Color online) Experimental data of a quantum dot 
with R = 70 nm. This figure is taken from d. The single- 
particle energy of nine consecutive states (labeled with red 
triangles and blue circles, respectivelt) is shown. The charac- 
teristic slopes of the dashed lines are about ±2.5meV/T. 

Employing the constant-interaction model, the ground 
state energy of an A'^-particle quantum dot can be writ- 
ten as 

E^.{N)^Y.^,{B) + ^~eNaV,. (12) 

i—l 

The single-particle energy ei{B) of the z-th particle is 
given by E{n,m,T). The second summand is the elec- 
trostatic contribution to the energy with the total capac- 
itance Cs of the experimental device which we assume 
to be independent of B. The ground state energy can be 
tuned by a gate voltage Vg. The lever arm a is deduced 
from Coulomb diamond measurements. The experiment 
described in 0] was done in the zero-bias regime; in 
other words we measured the chemical potential /iat = 
Egs{N) - Egs{N - 1) of the A^-th Coulomb resonance. 
Hence, the single-particle energy £n[B) can experimen- 
tally be determined by eN{B) — eaVg°^{N, B) + const. 
In Fig. [21 the constant part is subtracted so that con- 
secutive Coulomb resonances (labeled with red triangles 
and blue circles, respectively) touch each other in one 
point. Since there is no well-defined zero energy (the 
electron-hole transition could not be determined in the 
experiment), an arbitrary offset was subtracted. 

We compare the experimental data to our theoretical 
model. For this purpose, we zoom into a particular re- 
gion of Fig. [T] Since we are now interested in states with 
a constant number of particles, these are - in analogy to 
Fig. [2] - shown consecutively with red dashed and blue 
solid lines in Fig. [S] The kinks in the energy spectrum 
for a constant number N of particles occur whenever the 
A^-th particle changes its quantum state to stay in the 
lowest possible single-particle state. In Fig. [31 the slopes 
are constant and vary between 7meV/T and 12meV/T. 
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structures. However, apart from the reasonable agree- 
ment between theory and experiment, instructive mod- 
els like the one presented here can be solved analytically 
with the infinite-mass boundary and give an intuitive un- 
derstanding of the physics behind such systems. 

Acknowledgement: Financial support by ETH Ziirich 
and the Swiss Science Foundation is gratefully acknowl- 
edged. 

Note added: During preparation of this manuscript, we 
became aware of similar and more extensive calculations 
by P. Recher et al. 
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FIG. 3: (Color online) A close-up of the single-particle spec- 
trum shown in [1] For comparison to the experimental data, 
states with a fixed number of particles are consecutively de- 
noted by red, dashed and blue, solid lines. An arbitrary offset 
is subtracted as described in the text. Slopes marked by the 
two black, dashed lines vary between 7 meV/T and 12 meV/T. 



This is in reasonable agreement with our experimental 
data where we measured slopes of 2.5meV/T. The fol- 
lowing three reasons might limit the quantitative agree- 
ment between experiment and theory, i) The circular 
shape is a simplification of the experimental device, ii) 
Disorder is expected to play a significant role in graphene 
nanostructures. This is not included in the theoretical 
model, iii) Since the experiments were not performed in 
the single-electron regime, interactions effects should be 
included in a thorough theoretical analysis. Numerical 
tight-binding calculations or quasi-classical simulations 
are capable of implementing these aspects . 

III. SUMMARY 

In this paper, we derived an analytic expression for 
the energy spectrum of a circular graphene quantum 
dot which is exposed to a perpendicular magnetic field. 
The boundary condition we employed is the infinite-mass 
boundary introduced in [llj . This straightforward model 
is in good qualitative agreement with our recent experi- 
ments [sj. We discuss possible limitations of the model. 

The validity of the infinite-mass boundary condition 
for graphene is discussed in recent papers 15]. One needs 
to evaluate the consequences of different boundary con- 
ditions for energy spectra of confined graphene quantum 
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